Modeling the underlying mechanisms for organic memory devices: Tunneling, electron 

emission and oxygen adsorbing 
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We present a combined experimental and theoretical study to get insight into both memory and 
negative differential resistance (NDR) effect in organic memory devices. The theoretical model we 
propose is simply a one-dimensional metallic island array embedding within two electrodes. We 
use scattering operator method to evaluate the tunneling current among the electrode and islands 
to establish the basic bistable I-V curves for several devices. The theoretical results match the 
experiments very well, and both memory and NDR effect could be understood comprehensively. The 
experimental correspondence, say, the experiment of changing the pressure of oxygen, is addressed 
as well. 

PACS numbers: 73.23.-b, 72.80.-r, 73.40.Gk 



Memory devices based on metal-insulator-metal resis- 
tive structure with organic materials are promising in 
application^ The underlying mechanism in these devices 
is not very clear up to now due to the complicated phe- 
nomena observed, say for example, the memory effect and 
the negative differential resistance (NDR) in ON state, 
making an obstacle for controllable improvement of the 
performance Two different but correlated mechanisms 
have been addressed. The first one is the Simmons- 
Verderber (SV) model based on the space charge accumu- 
lated on some penetrated metallic islands^ This model 
could be proposed to explain NDR, which is also closely 
related to the memory effect j£ In addition to it, electron 
tunneling and Coulomb blockade were also proposed, ? 
indicating that the different spatial distribution of the 
space charge gave rise to the transition between ON and 
OFF state. The admittance spectroscopy provided an- 
other evidence to the SV models Whereas, the basic 
drawback of this mechanism is that, the OFF state is 
always observed, in practice, during the first scanning 
of bias voltage, while SV model predicts the sample is 
initially in ON stated The second mechanism is the fila- 
mentary conduction, stating that there would be some 
penetration of metal atoms into organic layer during 
the deposition of top electrode, which forms some con- 
ducting filaments under relative large bias voltage, and 
the trap-controlled tunneling between different filaments 
dominates the NDR effect The direct experimental ev- 
idence for this mechanism is from the transmission elec- 
tron microscopy image, which shows the metallic fila- 
ments are formed in ON stated Meanwhile, the tempera- 
ture insensitivity of ON state current^ the time depen- 
dent switch-on behavior and the scanning voltage 
and conductance dependence^ all go together to sup- 
port the filamentary mechanism. However, it is still hard 
to imagine the destruction and recovery of an atomic 
filament could respond so fast to the scanning voltage. 
Combining the above two mechanisms, we can get a com- 
mon idea that, during the deposition of the electrode, the 
metallic atoms could go deeply into the organic layer and 
form some islands in the film,— ^ which is crucial to the 



NDR effect no matter it originates from space charge or 
filaments. 

Another experiments were focusing on the role of in- 
terface. Some hot spots have been found on the surface 
of the device, and these spots exist only in ON state but 
disappear in OFF stated This in some sense implies the 
memory effect is most possibly caused by the build-up 
of some interfacial metallic tips, i.e., the OFF state from 
a flat metal-organic interface, while the ON state from 
some tips formed by penetrating metallic atomsJ^ The- 
oretical estimation shows that, at these tips the electric 
field could be as large as 5 x 10 6 V/cm, under which the 
field emission of electrons works. In order to verify this 
statement, we reported very recently an experiment on 
electric field and temperature dependence of switching 
time of the devices from OFF to ON stated The con- 
clusion there figured out the very possibility that metallic 
tips at the interface matter in the memory effect. 

Combining the above two considerations, it seems ap- 
propriate to say that, the NDR effect is from the metallic 
islands formed in the organic layer, while memory effect 
is from the electron emission at the interfacial metallic 
tips. Therefore, on the present stage, an in-depth and 
quantitative theoretical work responding to the relative 
experimental understanding is highly demanded. In this 
Letter, we will first establish a theoretical model to sim- 
ulate the metallic islands and tips in the organic memory 
devices and then calculate the tunneling current among 
the electrode and islands to study the basic characteris- 
tic of different devices. Then the experiment of changing 
the pressure of oxygen will be addressed to verify the 
corresponding theory. 

The model we consider in this work is briefly shown in 
Fig. [TJ The internal structure of the device is that, be- 
tween the two electrodes, there are some isolated metallic 
islands surrounded by oxygen molecules. The islands are 
formed during the deposition of top electrode, and the 
oxygen molecules are present when we do the measure- 
ment of pressure dependence, as will be discussed be- 
low. Because the islands are separated with electrodes, 
the electrons injected from the electrode to the islands 
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FIG. 1: Schematic of the structure of the system. The upper 
part is for the device structure and the lower one is for the 
energy diagram. See in the text for the detailed description. 



must first tunnel the potential barrier Uq formed by injec- 
tion barrier of normal organic molecules, which is around 
leV depending on the LUMO or HOMO of organics and 
Fermi level of metals. The oxygen molecules adsorbed 
by the metallic islands in the organic film will act as 
deep traps to capture electrons, such that the barrier 
should be largely increased due to the repulsion between 
the conducting electrons and trapped electrons.— Here, 
for simplicity, we have assumed that the distribution of 
adsorbed oxygen is an exponential decay from the top 
electrode with the decay length L Q , and the density of 
trapped electron is linear in space. The transfer integral 
J between neighboring islands is dependent on the dis- 
tance between these islands. Within this structure, the 
metallic islands could form some conducting channels for 
electrons, so that electrons could transport from one elec- 
trode to the other via quantum tunnelings. The trapped 
electron in adsorbed oxygen molecules will influence the 
potential profile of the islands, and when the bias voltage 
is strengthened, the effect of NDR should be founds On 
the other hand, when the device is switched to ON state, 
some metallic tips at metal-organic interface will form as 
discussed above, and then the injection barrier should be 
decreased due to the electron emission. In this sense, the 
current becomes larger and bistable effect emerges. 

Based on the above description, we write down the 
model Hamiltonian of the system as 



The first term for the electrode is 



H E - ^2(ta,k ± ^)Ci, fe Ca,fc, 



(2) 



where a = L/R is for left/right electrode, C\ k (C a ^) 
the electron creation (annihilation) operator of the k-th 
continuum state in the electrode, e a ,k the corresponding 
potential energy of that state, V& the bias voltage, and ± 
taking + for left and — for right. The second term of (JTJ) 
for the metallic islands is 



(3) 



where C\ (Ci) creates (annihilates) an electron at i-th 
metallic island, hi — CjCi, J the transfer integral be- 
tween islands, and n°{= fiNi) the charges adsorbed by 
oxygen molecules, which plays the essential role in this 
work. Wherein, iV,[~ exp(— i/L a )] the number of oxygen 
molecule adhering to the metallic islands, fa the filling 
factor, which will be assumed to be linear to i. The inter- 
facial injection potential induced by oxygen molecule is 
Uq, and the length of the system is L, which is not crucial 
when it is changed and will be set to 50 throughout this 
work. In all, there are three adjustable parameters in this 
theory: Uq is the interfacial injection potential, which de- 
termines the transition between OFF and ON state; J is 
the coupling between neighboring islands, which is deter- 
mined by how many islands are formed in the film, and 
we have neglected the disorder of J; L a is the penetra- 
tion length of oxygen molecule, which is determined by 
the concentration of oxygens and the trapped electrons 
in them as well. All these parameters could be directly 
determined by the experiments, as is shown later. 

The last term of (JTJ) is the coupling between electrodes 
and islands, which has the form 



Hc=J2 J o(Ct ik C 1/L + h.c.), 



(4) 
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where Jo is the transfer integral between electrode and 
island. This parameter is directly related to the spectral 
broadening induced by electrodes, such that one can use 
it as a scaling to all results and then compare them with 
the experimental measurements. 

Under a given bias voltage, a nonequilibrium steady 
state (NESS) could be achieved, and one can then calcu- 
late the corresponding I-V curve. To evaluate the tun- 
neling current of NESS, we utilize the scattering oper- 
ator method, which has been well developed by J. E. 
Han for the system of quantum dots*il The basic idea of 
the method is to derive the scattering operators from the 
Lippmann-Schwinger equation and then evaluate the cor- 
responding tunneling current i 1 ' 19 Basically, the method 
we use in this work is straightforward from Han's treat- 
ment, such that we would not go further to discuss the 
details of the method. 
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Based on the our model ([TJ , we first try to rebuild the 
basic experimental I- V curves of both ON and OFF states 
for different device structures. The details of experi- 
ments are as follows. All the materials used were com- 
mercially available without any further purification. The 
devices were fabricated on ITO-coated glass substrates 
with sheet resistance about 17fi/sq. Each substrate was 
heated to f50°C in air after solvent cleaning and then 
loaded into a high vacuum chamber (~ 2 x I0 _6 Pa). All 
the metallic and organic layers were formed by vacuum 
depositions. The deposition rate of Alq 3 , Mg and Al was 
about 0.8A/s, lA/s and I A/s respectively. The area of 
the organic memory device was 9mm 2 , defined by 3-mm- 
wide ITO stripes overcrossed by 3-mm-wide top metal 
stripes. The current-voltage (I-V) characteristics of the 
devices were measured by a Keithley 236 unit with the 
top electrode grounded. The measurements of electric 
properties of the devices were carried out in a glove box 
(O2 <fppm and H20<fppm) or a dewar (Janis Research 
Co. Inc. VTF-I00) under 10 _1 Pa pumped by a mechan- 
ical pump. 

In the left panels of Fig. [21 we show experimental 
measurements for three different structures by changing 
metallic electrodes. One can find some typical character- 
istics for these curves. For example, as the character of 
the memory effect, the OFF state current has a turn-on 
point at around 3V depending on the specific structure, 
while the ON state current goes up very quickly from 
OV. Both ON and OFF state current have the similar line 
shape, that is, a wavepacket-like shape in a certain region 
of voltage. Especially, the downward sides of both states, 
say the NDR region, overlap heavily. The ratio between 
ON and OFF state and also the absolute value of them 
depend on the specific structure. In the right panels of 
Fig. [21 we show the corresponding theoretical results. It 
is shown that, all the results are in good agreement with 
the experiments. By changing Uq, we have obtained the 
typical ON and OFF state current, respectively. That 
means, the memory effect is shown to correspond to the 
intcrfacial potential reduction, which is consistent with 
our previous experimental results^ On the other hand, 
the NDR effect has also been rebuilt by the theory, which 
is of importance in this type of memory device as we have 
mentioned. 

To explain the NDR effect we found, we summarize 
the whole physical process inside our theory as follows: 
Initially, when the scanning bias voltage starts from zero, 
all the energy potentials on each metallic island are the 
same, and no current is present. Then the voltage in- 
creases gradually, and the current follows to increase in- 
stantaneously. At this moment, the potential on each 
island shifts and starts to differ from each other. As we 
know, the largest tunneling current emerges once the po- 
tential of each island coincides. That means, under a 
certain voltage, depending on the potential drop among 
the islands, the current will meet the maximum value and 
then go down, so that NDR effect occurs. 

It is worth noting that, for each curve of the theoret- 
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FIG. 2: Experimental (left panels) and corresponding theo- 
retical (right panels) results of both ON and OFF I-V curves 
for different device structures: (a) ITO/Alq 3 (120nm)/Al; (b) 
ITO/Alq 3 (120nm)/LiF/Al; (c) ITO/Alq 3 (120nm)/Mg/Al; 
and different parameters: (d) J = 0.35eV, L — 0.1L; (e) 
J = 0.4eV, L = 0.LL; (f) J = 0.5eV, L a = 0.1L. 



ical results, the Uq change is about 0.7V from OFF to 
ON. Thus, if the metallic tip formed at the interface is 
estimated to be fOnm, the electric field it induces will 
be 7 x I0 5 V/cm, much smaller than the estimation from 
experimenters^ The condition for the bistable effect be- 
comes much softer in our theory. Meanwhile, the value 
of Uq for each device at OFF state is comparable to the 
injection barrier of electron, which is consistent with the 
fact that the metallic islands, as donors, matter in the 
whole process. 

To verify the present theory, we have to test it for in- 
dependent experiments. Here, we show the results for 
different oxygen pressures, in Fig. [3l The left panels 
are for the ON state current of different structures under 
various pressures. It shows that, with increasing pres- 
sure, the current reduces very quickly. As the OFF state 
current is also much smaller than the ON state, this re- 
duction of current is quite different from ON-OFF transi- 
tion. The reason is that, compared with Fig. [21 there's no 
turn-on point as the OFF state current shows, and also 
the downward sides separate largely. This indicates the 
oxygen plays a different role with the interfacial metallic 
tip, which matters in the ON-OFF transition. Another 
important point is that, we have also found (not shown 
here) the reduction of current is reversible, namely, when 
the pressure of oxygen is reset to a small value, the ON 
state current revives. More details of the experiment will 
be reported elsewhere. From this experiment, we can see 
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FIG. 3: Experimental (left panels) and corresponding theo- 
retical (right panels) results of ON state current dependence 
on the oxygen pressure and L for different device struc- 
tures: (a) ITO/Alqs (120nm)/Au; (b) ITO/Alq 3 (120nm)/Al; 
and different parameters: (c) J = 0.55eV, Uo = 1.2eV; (d) 
J = 0.35eV, Uo ~ 1.4eV. L ~ infinite means the oxygen 
molecule uniformly distributes in the organic film. 

that, the present of oxygen environment suppresses the 
memory effect heavily, which implies the essential role of 
oxygen in this type of memory devices. 

In a common sense, the oxygen penetrates into the or- 
ganic layer through the top electrode, due to the loose 
structure of organic materials. Hence, the penetration 
length L Q is intuitively proportional to the pressure of 
oxygen. To study this effect more comprehensively, we 
could then change the L a to check if the theoretical result 
behaves the same as experiment. As the right panels of 
Fig. [3] show, we have rebuilt the experimental data, ex- 
cept within the high bias voltage region, in which the IV 
curves of devices undergo a shift due to some unknown 
cause. Typically, both the turn-on point and the down- 
ward side of the curves are the same with those from 
experiment. The explanation of this finding is as fol- 
lows. When we change Uo, it means the potential drop 
among the islands becomes quicker, such that the maxi- 
mum value shifts to higher bias voltage. However, when 



we change L D , it means the potential drop in the film be- 
comes slower, but the injection at the bottom electrode 
becomes harder. Within this understanding, the maxi- 
mum point does not change very much, but the whole 
current curve goes down. This description could also ex- 
plain why in those materials with good conducting per- 
formance, such as pentacene and CuPc, the memory ef- 
fect is less robust That is, in these materials, the den- 
sity of electron trapped by oxygen is larger than normal 
cases. 

Finally, we would like to discuss the other type of or- 
ganic memory device, say for example, that based on 
the electrochemical reaction.- 5 It was proposed that, a 
creation-and-rupture process of redox-controlled metallic 
bridges is respondence to the memory effect. Although 
this electrochemical process is quite different from the 
formation of metallic tips, to simulate it one can still act 
to change the interfacial potential as we do in this work. 
It implies that, the microscopic origin for different types 
of organic memory devices might be different, but the 
underlying mechanisms are intrinsically the same. 

In summary, we have presented a theoretical model 
combining the effect of tunneling among metallic islands 
and interfacial electron emission due to the interfacial 
metallic tips to comprehensively investigate the under- 
lying mechanism of organic memory devices. We calcu- 
late several sets of theoretical results, which compare well 
with corresponding experimental measurements in differ- 
ent devices. Based upon this comparison between theory 
and experiment, we conclude that the two essential effects 
in these devices, say memory effect and negative differ- 
ential resistance, originate from the formation of metallic 
tips at interface and tunneling among metallic islands in 
the bulk, respectively. To further prove this conclusion, 
we show both theoretical and experimental results on the 
influence of oxygen, which shows a clear oxygen pressure 
dependence. Finally, we give an expectation on the in- 
fluence by changing the thickness of organic film. 
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